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Abstract 



We report the results of a molecular dynamics simulation of a supercooled 
binary Lennard-Jones mixture. By plotting the self intermediate scattering 
functions vs. rescaled time, we find a master curve in the /3-relaxation regime. 
■^ I This master curve can be fitted well by a power-law for almost three decades 

- ^ ' in rescaled time and the scaling time, or relaxation time, has a power-law 

C^ ' dependence on temperature. Thus the predictions of mode-coupling-theory 

H . on the existence of a von Schweidler law are found to hold for this system; 

— L , moreover, the exponents in these two power-laws are very close to satisfying 

Ph ', the exponent relationship predicted by the mode-coupling-theory. At low 

j^ I temperatures, the diffusion constants also show a power-law behavior with 

the same critical temperature. However, the exponent for diffusion differs 
from that of the relaxation time, a result that is in disagreement with the 
theory. 

PACS numbers: 61.20.Lc, 61.20.Ja, 02.70.Ns, 64.70.Pf 
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In the last few years a remarkable amount of activity has taken place in the field of the 
glass transition and the dynamics of supercooled liquids. This activity was mainly spawned 
by the development of sophisticated theories and novel experimental techniques with which 
the predictions of these theories could be tested. The most outstanding example of these 
theories is the so-called mode-coupling-theory (MCT) which was proposed by Bengtzelius, 
Gotze and Sjolander and, independently, by Leutheusser |l|]. The central point of MCT is the 
prediction of the existence of a dynamic singularity when the temperature of a supercooled 
liquid is lowered below a certain critical temperature T^. This point of view is in contrast to 
other theories that try to make a connection between the sluggish dynamics of glass forming 
materials in the vicinity of the glass transition and an underlying thermodynamic phase 
transition. 

For temperatures close to T^., MCT makes detailed predictions about the behavior of all 
correlation functions < X{Q)Y{t) > whose dynamical variables X and Y have a nonvanishing 
overlap with the density fluctuation pq for any wave vector q. These predictions have 
been verified in several experiments and computer simulations of supercooled liquids on a 
qualitative basis and, more recently, also on a quantitative basis. However, there are also 
experiments and computer simulations which show that the predictions of the theory are 
not always fulfilled 0,^. The reader can find a good introduction to the theory in the 
review articles by Gotze and Gotze and Sjogren |^ and most of the relevant references to 
the experiments and simulations in references [§,H. 

Since experiments and computer simulations have given mixed results as far as the valid- 
ity of MCT is concerned, more tests are clearly necessary. Moreover, quantitative tests of the 
theory would be useful in order to understand whether MCT merely suggests a method for 
performing scaling analyses of experimental data or whether it is also a correct and accurate 
theory of the exponents and exponent relations. Since MCT was originally developed for 
the description of simple supercooled liquids we have decided to test in a quantitative way 
the validity and applicability of the prediction of MCT for such a system and this letter 
reports some of our findings. 

Some of the most important predictions of MCT deal with the behavior of the correlation 
functions in the so-called /5-relaxation regime. At low temperatures the time dependence of 
correlation functions of supercooled liquids show a broad shoulder or even a plateau when 
plotted versus the logarithm of time. The approach to and the subsequent departure from 
this plateau defines the /5- relaxation regime [^. One of the main predictions of MCT for 
the behavior of the correlation functions in this /3-relaxation regime is the existence of a von 
Schweidler law. This law states that a correlation function 0(i) can be written in the form 

cP{t) = f-A{t/r{T)f . (1) 

Here / is the height of the plateau (the offset), often also called the nonergodicity parameter, 
and A and b are positive constants. A and / are predicted to be smooth functions of 
temperature if T > Tc and to be dependent on the type of correlation function investigated 
but the exponent b should be the same for all correlation functions. Furthermore the theory 
predicts that the relaxation-time t{T) shows a power-law divergence at T = Tc with an 
exponent 7 which is related to the value of b. Although the predictions of MCT concerning 
the von Schweidler law have experimentally been found to be true for some systems it has so 
far not convincingly been observed in computer simulations of simple liquids. However, for a 



lattice gas, such a behavior has been observed 0. But since a simulation of a different lattice 
gas has shown remarkable disagreement with the predictions of the theory 0], the application 
of the results of MCT to these kind of systems may be problematic. Therefore we decided to 
perform a large scale computer simulation of a simple liquid in order to investigate, among 
other things, the dynamics of supercooled liquids in the /3-relaxation regime. 

Molecular dynamics computer simulations are particularly well suited for testing the 
predictions of MCT since they allow the measurement of many different types of correlation 
functions. This freedom allows the performance of more stringent tests of the theory than 
would be possible for experiments. In addition, in computer simulations the measurement 
of the correlation functions is done in a very direct way. This is in contrast to experiments, 
where in most cases a considerable amount of interpretation and manipulation of the raw 
data has to be done in order to obtain the desired quantities. The main problems with 
molecular dynamics simulations are the limited size of the systems one can study and the 
time span that a simulation can cover. Since MCT is a theory of equilibrium it is of utmost 
importance to make sure that the simulated system is equilibrated at all temperatures 
investigated. A very recent simulation of a glassy system has shown that failure to equilibrate 
the system leads to a completely different relaxation behavior [0]. In order to overcome this 
problem we have performed a simulation that covers a time range that is more than an order 
of magnitude in time larger than previous molecular dynamics simulations that have studied 
relaxation in supercooled liquids. 

The system we are dealing with in this work is a binary mixture of classical particles. 
Both types of particles (A and B) have the same mass m and all particles interact by 
means of a Lennard- Jones potential, i.e. V{r) = 4e(((T/r)^^ — {a/r)^). The parameters e 
and a of the various interaction potentials were chosen as follows: caa = 1-0, ctaa = 1-0, 
Cab = 1-5, (Jab = 0.8, ess = 0.5, and cbb = 0.88. These potentials were truncated and 
shifted at a cutoff distance of 2.5cr. The number of A particles was 800, and the number 
of B particles was 200. These potentials are similar to the ones used by Stillinger and 
Weber for their simulation of amorphous Ni8oP2o [0- In the following, all quantities will 
be expressed in a system of units in which the unit of length is ctaa, the unit of time is 
(mcr^^/48e^^)^/^, and the unit of energy is eAA- The equations of motion were integrated 
with the velocity form of the Verlet algorithm with a step size of 0.01 and 0.02 at high and low 
temperatures respectively. This step size was small enough to restrict the fluctuations of the 
total energy below an acceptable level. The system was equilibrated at high temperatures 
and, by coupling it to a heat bath, subsequently slowly cooled down to low temperatures. At 
each temperature the system was equilibrated for a time which was equal to or longer than 
the time it took all the correlation functions investigated to decay to zero to within the noise. 
After equilibration, a long molecular dynamics calculation was performed for the calculation 
of the diffusion constant and the correlation functions. The length of this run was also equal 
to or longer than the time for the correlation functions to decay to zero, to within the noise 
of the calculation. The longest equilibrium runs were those for the lowest temperature and 
had a duration of 100,000 time units, which would correspond for an atomic liquid to a real 
time of about 10 ns. This is about an order of magnitude longer than previously reported 
comparable simulations. To improve the statistics further we determined the correlation 
functions for at least eight different starting positions in phase space. More details on the 
simulation will be presented elsewhere M. 



In the following we will mainly concentrate on the investigation of the dynamics in the 
/3-relaxation regime. Since the presence of activated, or hopping, processes will modify the 
prediction of the simple version of MCT (in which these effects are neglected) [^ one has to 
determine the importance of these processes for the system under investigation. By studying 
the self and distinct part of the van Hove correlation function we have found that on the 
time scale of the /3-relaxation there is no secondary peak in the self part and no peak at 
r = in the distinct part []TU[ . Thus we conclude that hopping processes are not relevant on 



the time scale of the /5-relaxation and can therefore be neglected in the analysis of the data 
1^. Thus the comparison of our results with the simple version of the theory is justified. 

The space Fourier transform of the self part of the van Hove correlation function gives 
the self intermediate scattering function Fs{q, t) ^^. Figure |l| shows the time dependence of 



Fs{q,t) for the AA correlation function for all temperatures investigated. The value of q is 
qmax, the location of the first maximum in the structure factor SA_A{q) for the AA correlation 
function. The following observations can be made: 1) For all temperatures investigated the 
correlation functions decay to zero to within the noise of the data. This is strong evidence 
that the runs were long enough to equilibrate the system. 2) For high temperatures (curves to 
the left) the correlation functions decay quickly to zero. When the temperature is lowered the 
formation of a shoulder at intermediate times (on a logarithmic time scale) can be observed. 
For the lowest temperatures (curves to the right) this shoulder forms almost a plateau and 
we can clearly observe the two step relaxation behavior observed for strongly supercooled 
liquids. 3) For times around 15 time units the correlation functions for low temperatures 
show a small bump. Lewis and Wahnstrom recently reported a similar observation in a 
computer simulation of the glass forming liquid ort/io-terphenyl [Q . These authors came 
to the conclusion that the location of this bump is a finite size effect. Therefore, and also 
since we have not observed this effect for the B particles, we won't discuss this feature here 
any longer. 

In order to test for the presence of a scaling behavior we plotted the data from Fig. |I| 
versus a rescaled time t/T{T). The value of the scaling time t{T) was chosen at each 
temperature such that Fs{q, r) = e~^. In Fig. ^ we show the resulting plot. We can clearly 
recognize the presence of a master curve. This master curve can be fitted very well by a 
power law of the form given by Eq. (|l|) with / = 0.783 and b = 0.488 ± 0.015. This fit 
is included in the plot as well. Since the fit is vahd for almost three orders of magnitude 
in rescaled time it is definitely significant. Thus we find a power law with an exponent b 
and an offset / that are independent of temperature. This is exactly the behavior predicted 
by MCT for the behavior of the correlation function in the later part of the /3-relaxation 
region, i.e. a von Schweidler law with a nonergodicity parameter / and an exponent b that 
are independent of temperature. 

In order to test whether the power-law observed in the /5-relaxation region is just the short 
time expansion of a Kohlrausch- Williams- Watt (KWW) law, 0(t) = /exp(— (t/r)^), which 
is known to be often an excellent approximation for the long time behavior of correlators 
in glassy materials, we fitted the long time behavior of our correlation functions with such 
a functional form. We found that such a fit is very good but that the extrapolation of the 
fit to intermediate times, i.e. to the /5-relaxation regime, is not good at all since the fit falls 
below the data. Thus we can conclude that the power-law observed here is not just the short 
time behavior of the KWW-law. 



MCT predicts that the scahng time r in Equ. ^scales with temperature as r oc (T—Tc)~'^. 
We fitted the relaxation time with such a power-law and found, for the A and B particles, 
for Tc the value 0.435. The exponent 7 is 2.5 and 2.6 for the A and B particles respectively. 
MCT also predicts a relationship between the exponents b and 7. The value found above 
for b, when combined with this relationship, predicts that 7 = 2.7, which is very close to 
what we actually found. Thus the exponent relationship of MCT is confirmed to within 
the precision with which we can obtain the exponents. In Fig. ^ we plot r~^ versus T — Tc 
(dotted curves). From this plot we recognize that the power-law behavior is observed for 
0.466 < T < 0.6. Although this range is not that large we recognize from Fig. ^ that it 
is only for this temperature range that the curves follow the master curve. Thus from the 
point of view of MCT we have a consistent picture with respect to this range. We also made 
a fit to the data with a Vogel-Fulcher law and found that this functional form gave a good fit 
to the data over a temperature range which is a bit larger than the one where the power-law 
is observed 0]. However, this observation is not in contradiction with MCT since the main 
point is that the theory works where it is supposed to work, namely close to T^. 

MCT predicts that the von Schweidler law should be present not only for q = qmax but 
all values of q and that the exponent b should be independent of q. We tested this prediction 
by computing the self intermediate scattering function Fs{q, t) for different values of q in the 
range 6.5o"^^ < 5' < 9.6cr^^, i.e. from q values a bit less than qmax, the location of the peak 
in the structure factor Saa^q), to values up to the first minimum of Saa{(i)- In Fig- i we 
plot these correlation functions as a function of t^ with b = 0.488. The part of the curves 
that are straight lines are power laws with exponent b. We recognize that this is the case 
for t^ lying between 2-3 and about 60. This corresponds to a time interval of 4-9 time units 
to 3600 time units. Thus we find that also this prediction of MCT holds for this system. 

We have also done similar calculations for the B particles and found a similar behavior 
as reported here for the A particles. The von Schweidler exponent b for the B particles 
is about 0.445 ± 0.015, and thus quite close to the one of the A particles. Since MCT 
predicts these two exponents to be the same this observation is also in accordance with the 
theory. Note that this kind of universality, i.e. the presence of a von Schweidler law for all 
types of correlators and that the von Schweidler exponent is independent of, or only weakly 
dependent upon, the correlator, is one of the main predictions of the theory. 

So far we have dealt only with the behavior of the correlators in the /3-relaxation regime. 
But MCT also proposes an intimate connection between this regime and the relaxation on 
the longest time scale, the so-called a-relaxation regime. For example, the theory predicts 
that the constant of diffusion D should show a power-law divergence at T^ with an exponent 
7 that is the same as the one for the relaxation times. We have determined D by fitting a 
straight line to the long time asymptote of the mean squared displacement. We have fitted 
D with a power-law of the form D oc [T — Tc)^ and for both types of particles found the 
critical temperature T^. to be 0.435. The exponent 7 was 2.0 for the A particles and 1.7 for 
the B particles. These results are shown in Fig. ^, where we plot D versus T — T^ in a double 
logarithmic way (dashed curves). Thus, in accordance with MCT, we find that the critical 
temperature for the constant of diffusion is the same as the one we found for the relaxation 
time T. However, the 7 exponents are slightly different for the diffusion constants of the 
two species, and both values are significantly different from the 7 obtained from the scaling 
time. These findings are in contradiction to the predictions of MCT. 



In the range were the data shows a power law behavior (i.e. T < 1.0) we also tried to 
fit it with a Vogel-Fulcher law, i.e. D ex exp{—E/{T — Tq)) and found that this type of fit 
is clearly inferior to the one with the power-law. 

In summary, the predictions of MCT concerning the von Schweidler law are fulfilled in 
quite an impressive way for the system investigated here, but the relationship predicted 
by MCT between the temperature dependence of the von Schweidler relaxation time and 
the temperature dependence of diffusion do not hold. For the /5-relaxation regime we have 
presented only the results of the simulation that deal with the departure of the correlation 
functions from the plateau. But to make a stringent test of the theory it is clearly necessary 
to check not only a few predictions of the theory but as many as possible. We have done 
this and the results will be presented in reference §] . These results seem to indicate f^ that 
the behavior of the correlation functions in the first part of the /3-relaxation, i.e. the critical 
decay, is not as well described by the theory as the second part which is described in this 
letter. Nevertheless, considering the approximations that have to be made in order to derive 
the statements of the theory, the accordance of its predictions with the results presented 
here is most remarkable. 
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FIGURES 

FIG. 1. Self intermediate scattering function Fs{q,t) for A particles vs. t. q = 7.251cj^^, the 
location of the first peak in the structure factor SAAiQ)- Temperatures (from left to right): 5.0, 
4.0, 3.0, 2.0, 1.0, 0.8, 0.6,0.55, 0.5, 0.475, 0.466. 

FIG. 2. Self intermediate scattering function for the A particles for q = 7.251(T^^ vs. rescaled 

time t/T{T) (solid lines). Dashed curve: Fit with von Schweidler law. Temperatures (from right 

to left) as in Fig. 2. 

FIG. 3. Inverse relaxation time (dashed curves, Tc = 0.430) and self-diffusion constant with 

error bars (dotted curves, Tc = 0.435) for A (circles) and B (squares) particles vs. T — Tc. Also 
shown are the fits with a power law (solid lines). 

FIG. 4. Self intermediate scattering function for A particles for 6.5cj^^ < g < 9.6a^^ (from 
top to bottom) vs. t with von Schweidler exponent b = 0.488. T = 0.466. 



